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Abstract. Gamma- Ray Bursts are good candidates of 
the "bottom up" scenario for the generation of the UHE 
Cosmic Rays. In the most discussed model of GRBs, 
namely the "fireball" model, a highly relativistic shock 
forms and seems capable of accelerating the cosmic rays 
up to the EeV range. However, only the first Fermi cycle 
produces a large energy gain to particles coming from the 
external medium. Thus, a complementary acceleration is 
proposed, downstream of the external shock, in the rela- 
tivistic plasma of the GRBs, where crossings of relativis- 
tic fronts are likely to occur. Both forward and backward 
fronts are necessary for the internal Fermi acceleration to 
work and the physical process that generates them is pre- 
sented. We found that there exists a relevant physical pro- 
cess similar to Brillouin backscattering that redistributes 
the incoming energy in the plasma shell. This redistribu- 
tion occurs through the generation of sound waves that 
heat the plasma shell and also through the generation of 
both forward and backward relativistic Alfven fronts that 
accelerate cosmic rays by the Fermi process. We show that 
this ensemble of processes is able to account for the gen- 
eration of UHE cosmic rays. 

There are two opportunities for these combined pro- 
cesses, first during the "primary" Gamma-Ray Burst 
where baryons are entrained by the relativistic pair wind, 
second during the predeceleration and deceleration stages 
when the fireball interacts with the interstellar medium. 
However, because of the synchrotron losses, only the sec- 
ond stage can produce the UHE cosmic rays. 

Key words: gamma-rays : bursts - cosmic rays - accel- 
eration of particles - plasma instabilities 



1. Introduction 

The origin of the cosmic ray population beyond lO^^ey is 
an important enigma of modern astrophysics which will 
hopefully be clarified by the Pierre Auger instrument. Be- 
yond this energy, the galactic magnetic field is unable to 



scatter cosmic rays and moreover supernovae remnants 
cannot accelerate protons. Pulsars could accelerate pro- 
tons up to IQ^'^eV. This population, that has an isotropic 
energy spectrum with a powerlaw index 3.1, probably 
comes from extragalactic sources. The "knee" of the spec- 
trum around lO^^eV is so smooth that the extragalactic 
contribution could be important at lower energy also. 

The photoproduction of pions by cosmic ray protons on 
the Cosmic Microwave Background, the GZK-effect, oc- 
curs beyond the energy threshold of 3 x lO^^eF and higher 
energy protons cannot come from sources located farther 
than lOOMpc (Aharonian & Cronin 1994 ). Few events 
beyond the GZK-threshold have recently been detected 
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by AGAS A (Ha yashida et al. |1998D and the "Fly's eyes" 
(Sokolsky 1998), and the Auger experiment will consider- 
ably improve the statistics of these events. These prelimi- 
nary results suggest that a new population of cosmic rays 
is observed at energies higher than 10^*ey, because the 
spectrum becomes harder. The origin of these few events, 
although some pointing ability of the instruments, remains 
very uncertain. The identified extragalactic sources within 
lOOMpc are very few. 

Vietri ( 1995| ) and Waxman (1995) showed that the rate 
of GRBs in the Universe {10~^yr~'^ Mpc~"^) and their en- 
ergy (10^^ — lO^^erg) make them very good candidates as 
sources of the UHECRs, if 10% of their energy is converted 
into UHECR energy. The most considered GR B mo del, 
namely the "fireball" model (Rees & Meszaros 1992), is 
based on a relativistic blast wave having a Lorentz fac- 
tor Fg larger than 10^, this high value being necessary for 
gamma photons to escape from pair production. A Fermi 
cycle at a relativistic shock can amplify the energy of cos- 
mic ray by a factor F^. Thus two Fermi cycles could be 
sufficient to reach the expected energy of the UHECRs. 
However, Gallant & Achterberg (1999) recently showed 
that only the first cycle can produce such amplification, 
the next cycles amplifying by a factor of 2 only. Moreover 
the escape probability is large (0.3 — 0.5). 

The GRB light curves indicate strong internal dis- 
turbances and they have been considered as the main 
cause of gamma-ray emission and a poss ible c ause of cos- 
mic ray generation (Waxman & Bahcall 1999| ; Daigne & 
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Mochkovitch 1995). In this paper we propose an inter- 
pretation of these disturbances in terms of relativistic hy- 
dromagnetic fronts, calculate their generation through the 
streaming instability caused by baryon loading, show that 
there is an efficient nonlinear generation of both forward 
and backward waves, estimate the amount of energy they 
extract from the fireball and calculate the acceleration of 
cosmic rays by their crossings. The baryon loading oc- 
curs first during the "primary" GRB when the relativistic 
pair wind entrains the debris of the progenitor, and then 
during the interaction of the fireball with the interstellar 
medium. Furthermore, the first loading is likely to occur 
in the "hypernova" scenario, however, this is less obvious 
in the "merging" scenario. 

2. Loading baryons 



The canonical fireball model (Rees & Meszaros 1992) 



presents a primary explosion during which the plasma is 
optically thick until {R^ ~ ct^), then an adiabatic and 
free expansion until the ambient medium starts to exert 
sufficient ram pressure; this defines a deceleration time td 
and a deceleration radius Rd (typically = 3 x 10~^td). 
The third stage is thus the deceleration one during which 
most of the energy of the fireball is dissipated in the inter- 
stellar medium, giving rise to the afterglow. We assume 
that the fireball is initially composed of e+ — pairs 
and that this ultra-relativistic wind of bulk Lorentz factor 
r, flowing along a strong poloidal magnetic field, is per- 
turbed by ambient baryons of mass density pi,. There are 
two stages of baryon loading, during the "primary" GRB 
and during the interaction with the interstellar medium. 



2.1. Hydrodynamics of the baryon loading stage 

The main features of the hydrodynamics of the relativistic 
fireball can be derived from the energy invariant as long as 
the radiation losses are negligible (assuming radial mag- 
netic field lines): 



+ 0'P)-i^d^'- 



(1) 



where e is the comoving energy-mass density, P the co- 
moving pressure, (3 the flow motion (in unit of light ve- 
locity), and 7 the Lorentz factor of the flow (7 = (1 — 
^2^-1/2^^ distributed around the typical bulk Lorentz fac- 
tor F. We assume that the plasma shell is composed of a 
low pressure baryonic component of mass M and a high 
pressure relativistic component containing pairs and cos- 
mic rays of average energy density in comoving frame 
e* = 3P. We approximate the energy by the following 
expression, assuming an appropriate self-similar evolution 
of the shell (which allows us to properly define the bulk 
Lorentz factor): 



E = MTc^ + 



(2) 



where Vb — AttR^SRq is the covolume of the shell, R{t) the 
shell radius, SRq its comoving thickness. During the ex- 
pansion, the energy in the shell is redistributed between ki- 
netic and internal energies through the competition of adi- 
abatic cooling and heating by the flux of incoming matter. 
Part of the mass flux increases the low pressure baryonic 
mass: M = (1 — a)pbFc47ri?^; the other part contributes to 
the high energy component: e* \heaUng= apbTc^A-KR"^ /Vq, 
assuming a heating length shorter than the fireball width. 
The adiabatic cooling is such that lad ~ 



j,diab- 



3 V(, 



Thus setting E = 0, after some simple manipulations, we 
obtain (with ao = 1 + a/3): 



aopb^TTR^c^r^ 
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It can easily be checked that, when there is no mass 
input, an adiabatic expansion occurs such that F = 
ro(Vb(i)/Vb(io))^/^ as long as MTc^ < E, followed by 
a free expansion F ~ E/Mc^ when the internal energy 
has become smaller than the kinetic energy. Now to an- 
alyze the effect of the mass sweeping, we assume a self- 
similar expansion such that dhiVo/dhit — x and define 
the dimensionless quantity •q{t) (within a coefficient of or- 
der unity) as the ratio of the fireball energy over the swept 
energy-mass of the ambient medium, measured in the co- 
moving frame, such that 



E^7^{t)pb--KR^T^^ . 
Then equation (||) leads to the differential equation 



(4) 



where r^o 
3 _ 
2 



dlni=2^(^°°-'^)' 

is the asymptotic value taken by ri{t), 
rjQ, the solution is 

770(^^0)'"- 



(5) 

namely 



f . With ??(io 



(6) 



The evolution starts in the adiabatic regime, 
Voit/toY^'^ and F cx t^^^. If rj reaches its asymptotic 
value, then the bulk Lorentz factor decreases according to 
the canonic law F cx t^^^^. However, the interest of this 



derivation is to stress that this Blandford & Mc Kee (1976) 
law does not mean that a free expansion has been reached 
and a strong relativistic shock set up. Indeed the energy of 
the fireball can still be dominated by the internal energy, 
so that a shock does not form, and the incoming baryons 
interact with the fireball plasma through a streaming in- 
stability. This is probably what happens during the pri- 
mary stage of baryon loading where a high entropy pair 
plasma entrains some baryonic mass, expected to be of the 
order 10~^Mq, coming from the debris of the progenitor. 
This interaction is probably more important than the final 
stage of interaction with the interstellar medium in gener- 
ating perturbations in the flow as described in subsection 
~3| and section |[ 
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Fig. 1. Evolution in time of the Lorentz factor F. The de- 
celeration giving rise to the afterglow starts at time <d, af- 
ter a free expansion involving the asymptotic bulk Lorentz 
factor Foo = EjMc^. 



Thus there are two stages of baryon loading (see 
Fig. 1^, the primary stage when the high entropy pair 
wind entrains a mass Mq ^ 10 "^M© coming from the de- 
bris of the massive star giving rise to the hypernova or 
of the merging compact objects, and a secondary stage 
corresponding to the sweeping of the ambient interstellar 
medium but that loads typically 10~^Mo only. 

2.2. The magnetic field and some scales 

We assume that this baryonic plasma flows along a strong 
poloidal magnetic field. This is the major assumption of 
this paper: B > Beq (we note Pp = [Beq/B]'^ < 1) and 
B > Bfn where B^q is the equipartition field correspond- 
ing to the plasma pressure in the fireball, and Bm is the 
equipartition field corresponding to the ram pressure ex- 
erted by the incoming baryons. 



B„ 



102( 



.1/2 r 

1 cm 10'* 



(7) 



Indeed, the poloidal magnetic pressure follows the same 
decrease as the relativistic pressure with distance; for con- 
ical expansion, Bp/Sn cx P cx i?^^. Moreover, a toroidal 
magnetic field is likely significant in the edge of the col- 
limated flow and can play a confining role. The decay of 
the toroidal field is weaker, because Bt cx R^^ in a coni- 
cal expansion. If the magnetic field dominates at the be- 
ginning, it remains so with the expansion. Moreover, a 
tiny initial toroidal component of the magnetic field can 
become dominant in the expansion and favors coUima- 
tion. Like in extragalactic jets, the poloidal field can be 
dominant along the axis of the flow and the toroidal field 
dominant in the edge of the flow. Thus we assume an in- 
ternal poloidal field such that Bp = Bd{R/ Rd)~'^ with 
Bd ^ Bm at Rd ^ lO^^cm and an external toroidal field 
Bt > Bd{R/Rd)^^; in particular at i?* (typically 10~^Rd), 
B ~ lO^G. Of course these estimates are sensitive to the 
model of magnetic field distribution, however the numbers 
we proposed are reasonable and help to explain out what 
happens. 



In the early stage the proton energy is severely lim- 
ited by synchrotron losses. Assuming an acceleration time 
tacc ~ i^i-L, where is the Larmor time, the maximum 
proton energy is such that the acceleration time is equal 
to the synchrotron time and gives 



= 2- 



101 



B 

Tg 



-1/2 



GeV. 



(8) 



Taking k ~ 10, the energy of lO^'^eV can be achieved 
only in the region where B is smaller than Bg ~ lO^G, 
which requires R > Rg — lO^^Rd- Bearing in mind that F 
reaches a value between lO^-lO^, the UHECRs (- 10'^°eV 
in the observer frame) are necessarily produced in the re- 
gion R> Rg. 

A neutrino radiation is produced during the early stage 
through pp-coUisions thanks to Fermi acceleration. In- 
deed, emax given by (||) exceeds IGeV when B < lO^^G, 
therefore where R > Rh ^ 10~^i?,. This low energy neu- 
trino emission ends when the pp mean free path becomes 
larger than 5R, thus when R = Rpp ~ 10^ Rh ^ R*. This 
stage of Fermi acceler ation supports the predictions made 
by Paczynski & Xu ( |1994 ), namely a neutrino emission 
between 30 GeV - TeV of global energy of a few percent 
of the fireball energy. 



2.3. Relaxation of the baryon stream in the fireball 

In this subsection, we examine how the backstream of 
baryons in the fireball, initially composed of e+ — pairs, 
undergoes a relaxation by triggering a beam instability in 
the pair plasma. 

The hydromagnetic perturbations so produced scatter 
the baryons that are then entrained and some of them ac- 
celerate to high energy. We presume that these perturba- 
tions are those revealed by the light curve and which also 
accelerate the electrons responsible for the synchrotron 
and inverse Compton gamma-emission. Indeed, initial per- 
turbations do not amplify in the expansion, they even de- 
cay (this differs from an expanding universe that is self- 
gravitating). Therefore we stress that these perturbations 
must be generated during the expansion by the appropri- 
ate instability. 

At the presumed energy of the particles, the Coulomb 
interactions are negligible (indeed, for relativistic elec- 
trons, the Coulomb cross section is as low as the Thomson 
cross section; as long as the shell is optically thin to Comp- 
ton scattering, the mean free path is larger than the shell 
width; this is of course more obvious when protons are 
involved in the collision) . Therefore the incoming baryons 
interact only with the magnetic field carried by the fire- 
ball. The magnetic field represents a more efficient obsta- 
cle if it is perpendicular to the flow, and the interaction 
starts in the form of an intense backward fast magneto- 
sonic wave. In a confined relativistic plasma, the fast mode 
propagates with a Lorentz factor = 7s7*, where the 
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sound Lorentz factor (corresponding to the parallel slow 
mode actually) 75 = a/3/2 and the Lorentz factor of the 
generalized Alfven waves (see Pelletier & Marcowith 1998) 
7, = (1 + :jg-)^/^. However, since we consider a collimated 
expansion, the most natural circumstance is the interac- 
tion along the parallel magnetic field within a wide solid 
angle around the axis. 

In this scope, the interaction scenario is the following. 
In the comoving frame, the proton back-stream, of velocity 
Vb — PbC and Lorentz factor Tt, along the field line gener- 
ates backward Alfven waves of vel ocity K = /3*c at a rate 
g given by (Marcowith et al. 1997 ; Pelletier & Marcowith 



1998 



9 = goTT VPb^-K 
1 6 



(9) 



where go = ^(2/3*) 



-1/2 



being the plasma frequency 



of the ions in the stream. This maximum growth occurs 
for the wave number 



7* 



Tbivb-V^) Tb c 



(10) 



When these backward waves have reached a high level, 
they scatter the baryons that are then entrained by the 
pair fiow. Thus, we presume that the large internal dis- 
turbances observed in the fight curve result from loading 
baryons in the fireball either at the early stage of the ex- 
plosion or during the "free" expansion stage. At the begin- 
ning of the fireball expansion in the interstellar medium, 
the diffusion length is not short enough compared to the 
narrow width of the shell. Then baryon entrainment starts 
when the level of the resonant Alfven waves (hereafter A- 
waves) is strong enough to get a diffusion length shorter 
than the shell thickness. Then the kinetic flux of the in- 
coming baryons is transformed into a backward fiux of in- 
tense A- waves. As long as these waves do not interact with 
the plasma shell (their wavelengthes are larger than the 
dissipation scale), the shell does not decelerate and just an 
electromagnetic wake is generated propagating towards to 
shell center but more slowly than the ultrarelativistic ex- 
pansion of the shell so that it seems to advance. Thus for 
a baryon mass loading rate Mb = AnR^pbTc, the energy 
flux, measured in the comoving frame, generated in the 
form of backward disturbances in the shell is: 



Mb 



(11) 



For a fireball of initial energy Eq, this baryon loading pro- 
duces only a weak shell deceleration as long as 



t < td 



MbTc^ 



(12) 



where is the deceleration time. In other words, the shell 
is still weakly perturbed as long as the mass loading rate 



is smaller than Mj, = pbTAirR^c where Rd = ctd is the de- 
celeration radius. There are in fact two deceleration radii 
corresponding to the two stages of baryon loading. 

The relaxation of the baryon stream in the fireball is 
described by the following simplified model coupling the 
relative level of Alfven energy density u with the relative 
motion of the stream /?{, — /3* with respect to the backward 
A- waves so-generated. The coordinate r denotes the radial 
distance from this external sheet, increasing r means ap- 
proaching the shell center and the "initial" condition is 

Tb{0) = r. 



du 
dr 
dTb 
dr 



2gu 

^lylSTb 



(13) 
(14) 



where in a quasi-linear approximation the slowing rate 
— vqu^ the more intense the waves the more efficient the 
scattering. The energy fiux conservation implies that 

MooP*pl b 

where Uoo is the asymptotic value of u, namely 
W./(BV87r) given by (|ll|) and thus Uoo = {B,n/B)^ < 1. 

The relaxation of the stream can be characterized by 
a single typical length Ir — K/5o- Fig- H illustrates the 
stream relaxation in the shell that occurs in few L. 



O.B - 



0.4 - 



Fig. 2. Stream relaxation and growth of the waves. The 
final relaxation stage, for /?</?*, is just a fast isotropi- 
sation involving plasma microphysics not detailed in this 
paper. The transition to the non-relastivistic regime gives 
rise to a stiff variation. 

Let us consider the interaction with the interstellar 
medium. The typical wavelength of the excited waves, 
measured in light seconds {I. sec), is such that 

B 




Ao ~ 10 



103'103G 

whereas the relaxation length is 



(16) 



(17) 
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This has to be compared with the shell width in the co- 
moving frame; at the deceleration time its estimate is 



(Rees & Meszaros 1992) 



=i 3 X 102( 



rib 



r 



1 cm ^ 



-5/3 



-2/3 



sec. 



10^^ erg' 

Let us consider now the interaction in the "primary" GRB. 
Expanding lO^^M© within a sphere of one light-second ra- 
dius (comoving) leads to a particle density of 10^^ cm~^. 
The plasma is now collisional, but the collective process 
previously developed provides a faster relaxation than 
both Coulomb and pp-coUision ones. The typical wave- 
length Aq — 10~"^"tht( iniiir^ )^"^ I. sec. and the relaxation 



length Ir ~ 10 



-10 r 
-9 r f 



rib 



Incidently, we 



mis " -a ) I. sec. 

10^° cm ^ 

note that although the time is quite short, the density is 
so high that the Lawson criterium for fusion is largely sat- 
isfied; moreover the thermal energy of the particles is high 
(e ~ 10^'^eV). Thus, alpha particles and neutrons should 
be produced and moreover, as mentioned previously, a sig- 
nificant neutrino emission. 

Once the relaxation is achieved, there are two ma- 
jor nonlinear interactions with the plasma shell through 
magnetosonic compression. First, intense backward A- 
waves that propagate almost parallel to the magnetic field 
produce fast magnetosonic compression governed by the 
Hada's system generalized to relativistic plasma by Pel- 
letier & Marcowith ( 1998] ). The transverse magnetic per- 
turbation b (reduced to the averaged field) exerts a pres- 
sure that produces fast parallel perturbed motion U|| (spe- 
cific momentum) that i s prop ortional to the parallel elec- 
tric field (see Pelletier 1999). This, parallel electric field 
is responsible for particle acceleration or stochastic heat- 
ing. For delocalized waves, particles that resonate with the 
parallel electric field produce the nonlinear Landau damp- 
ing of the A-waves. The so-produced E^^ efficiently injects 
electrons and positrons in the high energy population. 

There is also another strong nonlinear effect with 
the slow magnetosonic mode (S-mode) that efficiently 
backscatters the primary flux of A-waves. This is pre- 
sented in the next section. 

3. Brillouin backscattering 

As seen in the previous section, the flux of backward A- 
waves becomes stronger and stronger as long as the shell is 
still entraining baryons. If these waves are not be absorbed 
in the plasma shell, no deceleration nor heating of the shell 
would stem from the relaxation of the incoming baryons. 
Moreover, Fermi acceleration with these waves can work 
only if there are both forward and backward waves. One 
could expect some reflection of these waves in the internal 
edge of the shell; however, because they propagate at a 
velocity close to the light velocity, no significant change of 
impedance and thus no significant reflection occurs at the 



edge. Both absorption and backscattering efficiently de- 
velop by excitation of sound or slow magnetosonic modes. 

3.1. Parametric Brillouin instability 

A large fraction of this flux can be backscattered (thus 
producing a forward A-flux) by exciting the slow mode 
of the shell plasma. This is analogous to the Brillouin 
backscattering, but with A-waves instead of ordinary elec- 
tromagnetic waves. The analogous Raman scattering does 
not work with A-waves (three mode couplings does not 
work with A-waves). A mother A-wave (wq, ko) of relative 
amplitude 6o spontaneously gives rise to a slow magne- 
tosonic mode {uJs, fcs), but this S-mode is unable to carry 
the whole flux of energy-momentum and a backward (for- 
ward with respect of the shell expansion) A-wave is gener- 
ated with a lower energy-momentum {uj-, k^) such that: 



Wo 



T,d kn — fc. 



(18) 



For parallel propagating resonant waves, 



2fcn 



and k- = ko 



1 - Vs/V, 



one gets the most efficient rate of coherent wave decay 



(19) 



G decay = - ^CJg W_ ( ) 



(20) 



This result is a particular solution that can be derived 
from the complete analysis made by Champeaux et al. 
( 1999| ). For incoherent scattering, the decay rate is lower, 
proportional to 6§/2/3p instead of its square root. 

The S-modes are in turn absorbed by resonant interac- 
tion with "thermal" particles. This process is clearly the 
most efficient to transmit the momentum and energy of 
the baryon stream to the plasma shell. 

Now intense forward and backward flux of long A- 
waves are still remaining. There is no other way to ab- 
sorb them than cosmic ray acceleration through Fermi 
processes that involve resonant interaction between high 
energy particles and long A-waves. 



3.2. A toy model of the backscattering process 

In principle, solving the nonlinear system that governs 
the backscattering process allows us to predict how the 
dissipated energy is shared between heating the thermal 
particles of the shell and acceleration of cosmic rays. The 
detailed theory involves the numbers of quanta for each 
mode {N+^ N-,Ns) for the mother A-waves, the backscat- 
tered A-waves and the generated slow waves respectively) , 
a probability rate of scattering w and absorption rate 
(7+ , 7_ , 7s ) , all these triplets depending of the three wave 
vectors. In this paper we merely address a very simplified 
version of the system, similar to the system obtained with 
random phase approximation, but averaged over the spec- 
tral bands. We think it would be useful to solve such a toy 
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system to illustrate the process and to see whether reason- 
able estimates can be obtained by developing this method. 
For the sake of simplicity, we assume that backscattering 
occurs behind the stream relaxation sheet (in fact, it could 
even start in the relaxation sheet). The system is thus the 
following (in this notation a backscattered wave propa- 
gates outwards): 



d 

or 
d 

or 



-w{N+Ns + N+N_ - N_Ns) - 7+7V+(21) 
w{N+N, + N+N_ - N_N,) - 7_Ar_ (22) 
wiN+N, + N+N^ - N^Ns) - jsN, (23) 



The absorption rate 7s of the S-waves by the "thermal" 
plasma is always larger than the absorption rates 7+ and 
7_ of the A-waves by cosmic rays (gyro-synchrotron ab- 
sorption) ; moreover the S-waves have shorter wavelengths 
than the A-waves, the backscattered waves have even 
larger wavelengthes than their mother wave and thus even 
less absorbed. There are two scales in the system: a short 
characteristic scale associated to the growth of the S-waves 
I, = Vs/ijs - G), where G = w{N+{0) - iV_(0)) is the 
nonlinear scaterring rate, and a longer scale associated 
with the decay of the mother waves — T4(7s — G)/"fsG] 
Is ^ h when G <C 7s. This decay length U ~ fiX^/u. A 
ratio T of the incident flux is transmitted, a ratio R is 
reflected, a ratio A is absorbed and one has: 



T + R + A = l 



(24) 



and the absorption ratio is divided into a thermal one 
through S-waves, As, and a nonthermal one through A- 
waves, A*. Thus when K/7+ < <^^, T ~ 0. By solving the 
system, one obtains the backscattering rate: R — "-^-(°) 



One calculates the thermal absorption ratio: 

/o°° Is^sNsdr 



As 



V,LOoN+{Q) 



woAr+(0) 



(25) 



The absorption rate into cosmic rays is deduced by using 
(p^). Such a calculation is illustrated in Fig. ||. For reason- 
able absorption rates, we obtain R ~ 0.17, As ~ 0.53 and 
Atf ~ 0.30. Of course, a more realistic calculation should 
be made to predict these ratios, but this simple model in- 
dicates that, through the backscattering process, a sizable 
fraction of incoming energy can be converted into cosmic 
rays. In particular, a more detailed calculation would take 
into account that the longer wavelength modes are less 
damped and thus can be transmitted behind the shell. 

The combined process of stream relaxation and Bril- 
louin scattering followed by energy absorption by particles 
occurs first during the "primary" stage of the relativistic 
expansion when the pair wind entrains ambient baryons 
coming from the progenitor debris. As argued previously, 
this interaction is likely to occur without shock formation. 
Then after a new adiabatic and free expansion, a new stage 



g 0.2 




4 6 
distance 



Fig. 3. Solutions of the backscattering system. By order of 
decreasing line thickness we show 7V+/7V+(0), iV_/7V+(0) 
and A^s/-^^+(0). The distance unit is the scattering length 
defined by V^, /wN+{Q). In this example 7+ — 0.1,7- = 
0.05 and 7s ^ 0.69. 



of interaction occurs with the interstellar medium. The ra- 
tio of the energy dissipated by the cosmic rays acceleration 
over the initial fireball energy is E^r/E = A,,R{tY / R^ be- 
fore the deceleration time; thus, R = 0.3Rd, about 10% 
of the fireball energy in converted into cosmic rays. When 
the fireball radius reaches the deceleration radius, a strong 
collisionless shock is set up. However it can be considered 
as a hardening of the previous process. Indeed, incoming 
protons are reflected in the relaxation front giving rise of 
the external shock, whereas the backward magnetosonic 
wavefronts turn to internal shocks. The previous descrip- 
tion applies during the predeceleration stage. However, al- 
though the theory of beam-plasma instability is no longer 
relevant when the shock is set up, the Alfvcn fronts are 
still behind the shock and still accelerate cosmic rays effi- 
ciently. 

4. Fermi acceleration of UHE Cosmic Rays 

As stressed by Pelletier ( [1999 ), intense forward and back- 
ward A-waves of long wavelength (delocalized or localized) 
are necessary to accelerate UHE Cosmic Rays. The intense 
localized fronts propagate with a nonlinearly modified ve- 
locity such that 7, 1-^ 7„i = + 5j with Sj/j^, ~ 
{bm being the maximum reduced amplitude of the per- 
turbed magnetic field of the front measured in its comov- 
ing frame). These intense wave packets scatter particles 
of Larmor radius smaller than their width or wavelength 
in their frame in a few gyro-periods. In a scattering time, 
the particles energy gains a factor 7^ and the momenta are 
concentrated in the front cone of half-angle 1 /7* . Further 
acceleration requires other fronts propagating in the oppo- 
site direction, because wavefronts propagating in the same 
direction at almost the same velocity 14 tend to roughly 
isotropize the distribution of interacting particles with re- 
spect to their comoving frame. Thus it is crucial to get 
both forward and backward fronts to accelerate cosmic 
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rays and the Brillouin backscattering process of section 
^, that generates longwavelength perturbations in a time 
~ G^^, is the appropriate and efficient solution. The rel- 
ativistic regime of Fermi acceleration is deeply different 
from the nonrelativistic regime, not only because of the 
anisotropy effect but also because the acceleration time 
scale can become shorter than the scattering time scale. 
This is a crucial advantage that is illustrated by the fol- 



lowing result (Pelletier 1999) obtained in the case of a 



permanent flux of both forward and backward waves, the 
energy gain during a time At, smaller than the scattering 
time, is such that 



< Ap' >= p^Pl-tlusAt , 



(26) 



where Vs is the scattering frequency. Moreover, because 
B > Seg, 7* > 2, and the energy jump at each scattering 
is large and the statistic evolution cannot be treated by a 
Fokker-Planck description. 

Let us give some more details about the localized 
fronts. Intense longwaves tend to self-organize in for- 
ward and backward relativistic fronts of amplitude 5„j ^ 
(r*/^)~^/^, where ^ is the front width and r* = (< 7^ > 
/"f)Vi,/u!c is the minimum scale of relativistic MHD (Pel- 
letier & Marcowith 1998); it is worth noting that the 
perturbation amplitudes appear times larger in the 
observation frame. At each crossing of these fronts, the 
gain is by a factor 7^ within few gyro-periods. Such lo- 
calized fronts can be approximately described as solitons 



(Pelletier & Marcowith |l998t Pelletier |1999D , they dif- 
fer from ideal solitons because of the Fermi acceleration 
that produces a kind of Landau-synchrotron damping of 
them. When ideal solitons cross each other they do not 
destroy, whereas damping destroys them and their com- 
plete absorption corresponds to the strongest efficiency of 
the Fermi pro cess. A similar idea was used by Daigne & 
Mochkovitch ( |1998| ) with internal shocks to explain the 
gamma emission. Here we interpret the internal shocks as 
internal Alfven fronts which can contribute to short scale 
variations in the light curve, and possibly more easily than 
pure hydrodynamic shocks. 

For a bulk Lorentz factor of the shell of 10"^, particles 
in GRBs have to reach lO^^eF in the comoving frame in 
order to supply the UHE Cosmic Rays population beyond 
the GZK threshold. We consider the interaction with the 
interstellar medium only because, as previously showed, 
synchrotron losses prevent UHE cosmic ray production, 
inside a radius Rs (subsection |2.2| ). Assuming that the 
largest wavelengthes are one tenth of the effective shell 
width (limited by causal connection), a magnetic field just 
larger than the equipartition value at deceleration radius 
{Beq ^ IQ^G say) is enough to get Larmor radii of that size 
(i.e. ~ 10~^pc), which corresponds to an energy of order 
W^eV. With 7 ~ 10^, the scale 10 I. sec. Therefore, 
few fronts having a width of 1 — lOr* that cross each other 
within a shell width of 3 x IQ^l.sec. can produce protons 
of IQ^^eV within one second with respect of an observer 



frame during the interaction with the interstellar medium. 
When the strong relativistic shock has formed (i? > Rd), 
external particles gain energy by a first Fermi half-cycle 
by a factor F^. They are then injected in the relativistic 
A-fronts behind the shock and suffer further acceleration 
up to the ultimate energy. 

The local distribution is more likely quasi monoener- 
getic rather than a powerlaw. However the local charac- 
teristic energy is proportional to the product BR and thus 
the global distribution reflects the distribution of the prod- 
uct BR] as stated by Pelletier ( |1999[ ), a distribution close 
to e^^ is obtained for B ex with m close to 2. 

Our main conclusion is that Gamma-Ray Bursts are 
capable of producing the UHE Cosmic Rays through a 
Fermi process with relativistic Alfven waves. These in- 
tense waves are generated by the two stages of baryon 
entrainment, first during the primary GRB stage, second 
during the interaction with the interstellar medium. The 
Fermi process works because of the Brillouin backscatter- 
ing process that turns out to be appropriate and efficient, 
moreover it allows the heating and deceleration of the shell 
plasma by the incoming flux. Although the primary stage 
does not produce UHE cosmic rays because of synchrotron 
losses, the Fermi acceleration allows the maintenance of a 
high energy proton population between Rh and Rpp that 
can emit a significant low energy neutrino flux of the or- 
der 10~^E. A high energy neutrino flux can be generated 
through the p7-process by UHE cosmic rays at the end of 
the free expansion, typically a sizable fraction of the UHE 
cosmic ray energy in the form of neutrinos 



(Rachen & Meszaros 199S). 
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